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ABSTRACT 



When image data are processed to generate pseudo-tone 
patterns, a predetermined threshold is modified according to 
at least two periodic functions with different spatial periods. 
The picture elements are compared with the modified 
threshold, and the comparison results are used to generate an 
output image. The periodic threshold modifications cause 
dots to appear preferentially in regular patterns, thereby 
reducing dot noise, and providing improved control over the 
overlaying of dot patterns in different color planes in a color 
image. 

28 Claims, 16 Drawing Sheets 



31 -J 



WEIGHTING 
TABLE 



W1 



Tl 



T2 



T3 



T4 



rW2 



K3 



& 

M2 | 



W3 



& 



K4 



M3 



32 

jL 



EH J=pU 




-8A 

RING COUNTER 



12/04/2002, EAST Version: 1.03.0007 



U.S. Patent Oct 1, 2002 Sheet 1 of 16 



US 6,459,817 Bl 



FIG.1 



BLACK 
□ WHITE 



SI Imax ' 
2 0 



0 

S2 -Imax 
2 



Imax 



Imax/2 
0- 



EFl 



(Imax/2) +S1+S2 

/ 



1_ 



IMAGE 



12/04/2002, EAST Version: 1.03.0007 

\ 



1 t 



U.S. Patent Oct 1, 2002 Sheet 2 of 16 US 6,459,817 Bl 



FIG.2A 



■128. 
42 - 



2x2 



-43 



FIG.2B 



FIG.2C 



-128 


127 


42 


-43 



PM2 



4X4 



8X8 





















-] 








i; 




















































4 


2 




































-r 





PM3 

/ 



FIG. 2D 









































































































-1 


128 














15 


!7 














































































































































































































































4 


2 
















43 









































































































































PM4 



/ 



16x16 



12/04/2002, EAST Version: 1.03.0007 



U.S. Patent Oct 1,2002 Sheet 3 of 16 



US 6,459,817 Bl 




12/04/2002, EAST Version: 1.03.0007 



U.S. Patent Oct 1, 2002 Sheet 4 of 16 



US 6,459,817 Bl 



FIG.5A FIG.5B FIG.5C 



WHITE -r 255 WHITE -r 255 



Th' 
ATh 

Th = Imax/2 Th 



BLACK -L 0 



--128 



Th- 
Th'- 



BLACK -LO 



128 
ATh 



255-T + 127 



128- -0 



+42 



---43 



(H--128 



FIG. 6 



WEIGHT 



W4 W3 W2 




85 170 255 

GRAY LEVEL 



12/04/2002, EAST Version: 1.03.0007 




12/04/2002, EAST Version: 1.03.0007 



U.S. Patent Oct 1, 2002 Sheet 6 of 16 



US 6,459,817 



y < 



-128 



42 



FIG.8 

T1 



2x2 
T2 



127 
-43 



| - 128 1 


1 m 1 


1 42 1 


| —43 | 



4X4 
T3 













4- 






- 


L28 






















































4! 


I 
































T 





8X8 
T4 









































































































-] 
















12 


7 














































































































































































































































4! 


I 
















13 









































































































































16x16 



12/04/2002, EAST Version: 1.03.0007 



U.S. Patent 



Oct 1, 2002 Sheet 7 of 16 



US 6,459,817 Bl 



FIG.9 



x< 



1 

2 
3 
4 



1 

2 
3 
4 



FIG. 10 

PRIOR ART 



El 



Th 



E1- 



E2 N \fl 



E5 



E6 





A 

E3 
H2_ 



■E5 
E4 



2 



E7 



JET 



E7 



E6 



P1 P2 P3 P4 P5 P6 P7 P8 



□ 



12/04/2002, EAST Version: 1.03.0007 



U.S. Patent 



Oct 1, 2002 Sheet 8 of 16 



US 6,459,817 Bl 



FIG. 1 1 

PRIOR ART 




P1A 



FIG. 12 

PRIOR ART 



E (x - l,y - 1) A (x- l,y - 1) = 1/16 

E(x,y-1) 

A(x,y-1)=5/16 

E(x+l,y-l) 
A(x+l,y-l)=3/16 




E(x-l,y) 



P2(x, y) 
A(x-l,y)=7/16 



12/04/2002, EAST Version: 1.03.0007 



U.S. Patent 



Oct 1, 2002 Sheet 9 of 16 

FIG.13 



US 6,459,817 Bl 



Q START 



STEP 1 



STEP 2 



STEP 3 



INITIALIZE ERROR 
MEMORY 



CALCULATE E(x, y) 



STEP 4 



I'Gc y)^ 
Kx, y)+E(x, y) 



STEP 5 



GET Kl- 


-K4 


FROM 


TABLES 


Tl- 


-T4 



STEP 6 



CALCULATE Th' (x, y) 



STEP 7 




STEP 8 



OUTPUT SYMBOL 



STEP 9 



0 



OUTPUT SYMBOL 
"0* 



E(x, y)^ 

T (x, y)-Imax 



STEP 10 



E(x, y)^ 

I'(x, y)-Imin 



STEP 11 



12/04/2002, EAST Version: 1.03.0007 



U.S. Patent Oct 1, 2002 Sheet 10 of 16 US 6,459,817 Bl 



FIG. 14 



© 



STEP 12 




STEP 14 



YES 



STEP 13 



x— x+1 




YES 



STEP 15 



x— 0 

y^y + 1 



( END ) 



® 



12/04/2002, EAST Version: 1.03.0007 



U.S. Patent Oct 1, 2002 Sheet 11 of 16 US 6,459,817 Bl 




12/04/2002, EAST Version: 1.03.0007 



U.S. Patent Oct 1, 2002 Sheet 12 of 16 US 6,459,817 Bl 

FIG. 16 



PIXEL 



-n- 



m 



m 



A 



-n- 



-n- 



BLOCK 



FIG. 17 



Ab(xb-l,yb-l) 

Ab(xb,yb-1) 



B,LOCKS 











* - 


7^: 



Ab(xb+l,yb-l) 



Ab(xb-l,yb) 



12/04/2002, EAST Version: 1.03.0007 



U.S. Patent 



Oct 1, 2002 Sheet 13 of 16 

FIG.18 



US 6,459,817 Bl 



( START ) 



STEP 101 
STEP 102 

STEP 103 
STEP 104 

STEP 105 

STEP 106 
STEP 107 
STEP 108 



xb^O 
yb^O 



INITIALIZE 
ERROR MEMORY 








CALCl 
Eb(xb, 


JLATE 
yb) 






i— 1 








I' (xb, yb, i) 
— I(xb, yb, i) 
+Eb(xb, yb)/L 






GET K1-K4 FROM 
TABLES T1-T4 






CALCU1 
ATh(xl 


LATE 
b, yb, i) 






CALCU 
Th' (xb 


LATE 
, yb, i) 



XD, 3 



0 



® 



12/04/2002, EAST Version: 1.03.0007 



U.S. Patent Oct 1, 2002 Sheet 14 of 16 US 6,459,817 Bl 

FIG.19 



STEP 109 




STEP 112 



STEP 110 



OUTPUT SYMBOL "1" 



STEP 111 



STEP 114 



OUTPUT SYMBOL "0" 



E(xb, yb, i) 

— I' (xb, yb, i)-Imax 



STEP 113 



E(xb, yb, i) 

— -I'(xb, yb, i)-Imin 



i — i + 1 



STEP 115 




YES 



STEP 116 



STEP 117 



STEP 119 



Eb(xb, yb) 

— ZE(xb, yb, i) 



(b 




YES 



STEP 



xb— xb+1 



YES 



18 



( END ) 



~] STEP 



120 



xb — 0 
yb — yb+1 



12/04/2002, EAST Version: 1.03.0007 



U.S. Patent Oct 1, 2002 Sheet IS of 16 US 6,459,817 Bl 



FIG. 20 



IMAGE 



g3 

Th3 

S2 
Th2 

gl 
Thl 

g0 



OUTPUT 




12/04/2002, EAST Version: 1.03.0007 



U.S. Patent Oct 1, 2002 Sheet 16 of 16 US 6,459,817 Bl 




12/04/2002, EAST Version: 1.03.0007 



US 6,459,817 Bl 

1 2 

IMAGE-PROCESSING METHOD AND BRIEF DESCRIPTION OF THE DRAWINGS 

APPARATUS GENERATING PSEUDO-TONE , . . , . 

PATTERNS WITH IMPROVED REGULARITY 1D auacnca a^wings. 

FIG. 1 shows a one-dimensional example of the modifi- 

BACKGROUND OF THE INVENTION 5 cation of a threshold according to the present invention; 

The present invention relates to a method and apparatus FIGS. 2 A, 2B, 2C, and 2D show threshold pattern matri- 

for processing a gray-scale image or a color image to CCS cm P lo y cd in embodunents of the invention; 

generate pseudo-tone patterns with improved regularity. FIG - 3 illustrates a two-dimensional penodic function 

Pseudo-tone pattern generation is method of representing *V**** of ** P att ^ matrix in FIG. 2A; 

input gray levels or colors that cannot be directly expressed F1G * 4 illustrates a two-dimensional periodic function 

by an output device. Pseudo-tone pattern generation is comprising repeated copies of the pattern matrix in FIG. 2B; 

basically a type of quantization process, and is carried out by FIGS. 5 A, 5B, and 5C illustrate threshold modifications; 

comparing the image data with one or more threshold FIG 6 illustrates weighting functions used in the embodi- 

values. Loss of colors or gray levels can be avoided by 1S ments; 

diffusing the quantization error of each picture element or FIG ' 7 ^ a block ^ of a ^ cmbodimcnt of the 
pixel to some of the neighboring pixels. A more detailed invention- 
description of this error-diffusion process will be given later. *„ 4 . « c t1 ^ ^ 
r , • iL c • • * FIG. 8 illustrates a method of addressine tables Tl to T4 
For a gray-scale image, the result of error diffusion is that an . " ^ W1C * 

intermediate gray level is simulated by a pseudo-tone pattern 20 m 

of dots quantized to levels above and below the intermediate nG - 9 illustrates an alternative method of addressing 

level. The general term for this type of simulated gray-level tables Tl to T4 in FIG. 7; 

expression is dithering. Dithering can also be performed FIG. 10 illustrates conventional error diffusion in one 

directly, by replacing a single threshold with a matrix of dimension; 

differing thresholds. 25 FIGS. 11 and 12 illustrate two-dimensional error diftu- 

A typical use of error diffusion and dithering is to enable sion; 

a printer that can print only black and white dots to express FIGS. 13 and 14 are a flowchart illustrating the operation 

gray levels as patterns with appropriate densities of black 0 f the first embodiment; 

dots. A problem is that the dot patterns tend to exhibit FIG. 15 is a block diagram of a second embodiment of the 

random clumping and other random effects, generally 30 invention* 

referred to as dot noise, that do not look well when printed. ___ ' . . 

In many cases, a more regular dot pattern is desirable. For u nG ' " muS ™ 6S ^ dlVIS10n ° f m lnto blocks m 

1 j * »* *uu * the second embodiment; 

example, a pseudo-tone pattern expressing a fifty-percent ' 

gray level is preferably a checkerboard pattern of black and FIG ; 17 illustrates block error diffusion in the second 

white dots, rather than a random pattern. 35 embodiment; 

Error diffusion and dithering are also employed in color FIGS. 18 and 19 are a flowchart illustrating the operation 

printing. Here a further problem is that the random place- of me second embodiment; 

ment of the dots frustrates control over the way in which FIG. 20 illustrates thresholds employed in a third embodi- 

different color planes are overlaid. ment of the invention; and 

40 FIG. 21 is a block diagram of a fourth embodiment of the 



SUMMARY OF THE INVENTION 



invention. 



It is accordingly an object of the present invention to 

increase the regularity of the pseudo-tone patterns in a DETAILED DESCRIPTION OF THE \ 

quantized image. INVENTION 

Another object of the invention is to provide control over 45 Embodiments of the invention will be described with 

the relative positions of dot patterns in different color planes reference to the attached drawings, in which like parts are 

in a quantized color image. indicated by like reference characters. 

The invented method of processing image data comprises First? a simplified one^dimensional description of the 

the steps of: 5o modification of a threshold value by two periodic functions 

(a) modifying at least one predetermined threshold will be given. Consider a row of sixteen pixels in a gray- 
according to at least two periodic functions with dif- scale image having a uniform gray level equal to three- 
ferent periods; fourths the maximum level, the maximum level (1^ 

(b) comparing the image data with the modified threshold; representing white. When the gray-scale image is converted 
and 55 to a bi-level image, four of these sixteen pixels should be 

(c) generating output image data from the comparison made black, and the remaining twelve should be made 
results, the output image data being limited to a certain white. The basic threshold level (Th) is half the maximum 
number of pixel levels. level A first periodic function SI, which generates 

The invention also provides image-processing apparatus a periodic pattern, and a second periodic function S2, having 

for carrying out steps (a), (b), and (c). 60 ^ values illustrated in FIG. 1, are added to the basic 

The periodic modifications of the threshold cause dots to threshold level to obtain a modified threshold (I ma;e /2+Sl+ 

appear preferentially in predetermined positions, which have S2) that varies as shown at the bottom of FIG. 1. The image 

a regular, periodic pattern. Dot noise is thereby reduced. data are compared with the modified threshold, producing 

When a color image is processed, the alignment of the the desired four black pixels, as shown at the top of FIG. 1. 

regular patterns in different color planes can be controlled so 65 In two dimensions, a periodic function can be built up 

that the dot patterns are overlaid on one another in desirable from repetitions of a threshold offset pattern matrix, referred 

ways when, for example, the output image is printed. to below simply as a pattern matrix. FIGS. 2A, 2B, 2C, and 
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2D show four examples of pattern matrices with different The invention is not limited to pattern matrices of the type 

sizes. Each of these pattern matrices can be repeated to cover shown in FIGS. 2A, 2B, 2C, and 2D. For example, without 

a two-dimensional image plane, producing four periodic changing the algebraic signs of the values in the lower 

functions with different repeating periods. FIG. 3 shows the quadrants, these values (42, -43) can be replaced with the 

plane covered by copies of the pattern matrix PM1 in FIG. 5 values (127, -128) used in the upper quadrants, so that an 

2A, the solid lines indicating the matrix boundaries and the even positive-negative balance is also obtained in each 

dotted lines indicating cells within the matrices. The peri- vertical column. 

odic function is obtained by placing the values shown in The offset ATh is obtained separately for each pixel by 

FIG. 2A in the cells. FIG. 4 shows the plane covered by multiplying the values of the four periodic functions at the 

copies of the pattern matrix PM2 in FIG. 2B, giving a similar 10 pi^l position by four weighting functions Wl, W2, W3, 

periodic function. W4> ^ adding ^ resulting products. The weighting 

Each pattern matrix is divided into four quadrants, and has functions are functions of the gray level of the pixel. Each 

a constant value within each quadrant. The periodic function weighting function has non-zero weight values following a 

formed from the two-by-two pattern matrix PM1 thus has sine curve in a certain range of gray levels, and is zero 

the highest spatial frequency. The periodic function formed 15 outside this range. 

from the sixteen-by-sixteen pattern matrix PM4 has the Referring to FIG. 6, the weighting function Wl used with 
lowest spatial frequency. me highest-frequency periodic function (pattern matrix 
In the present invention, a desired periodic function is PM1) has zero weight for gray levels up to one hundred 
decomposed into several basic periodic functions in the seventy, and rises from zero to one in the range of gray levels 
same way as in a Fourier expansion. In the present example, 20 from one hundred seventy to two hundred fifty-five. The 
the basic functions are PM1, PM2, PM3, and PM4. The weighting function W2 used with the second-highest- 
basic functions are weighted and summed to build up the frequency periodic function (PM2) has positive weights for 
desired periodic function. The weighting coefficients depend gray levels from eighty-five to two hundred fifty-five, with 
on the input gray level. a peak weight (1) at one hundred seventy. The next weight- 
It will be assumed in the first embodiment below that the 25 ing function W3 is positive in the range of gray levels from 
gray scale varies from zero (black) to two hundred fifty-five zero to one hundred seventy, with a peak weight (1) at 
(white), and that image data having this gray scale are to be eighty-five. The weighting function W4 used with the 
converted to bi-level data. The basic conversion threshold is lowest-frequency periodic function (PM4) falls from one to 
one hundred twenty-eight (128, approximately equal to zero in the range of gray levels from zero to eighty-five, and 
255/2). This basic threshold (Th) will be modified by adding 30 is zero for gray levels above eighty-five. The positive weight 
a positive offset ATh as shown in FIG. 5 A, or a negative values of adjacent weighting functions thus overlap, with 
offset ATh as shown in FIG. 5B, to obtain a modified one weighting function increasing while the other weighting 
threshold Th'. The image data will be compared with the function decreases. The sum of all the weights at any given 
modified threshold. gray level is one, or a value close to one. (The sum of sine 
The modified threshold Th' is within the range from black 35 factions with a phase difference of ninety degrees varies 
to white if the offset ATh is within the range from minus one from one t0 approximately 1.4, which is close enough to one 
hundred twenty-eight to one hundred twenty-seven (-128 to for P reseDt purposes.) 

127). These two values are used in the two upper quadrants The process of multiplying the above periodic functions 

of each of the pattern matrices shown in FIGS. 2A, 2B, 2C, by weighting functions and summing the results is a type of 

and 2D. Values substantially equal to one -third of these Fourier integration process. 

values are used in the bottom two quadrants; these values Next, a first embodiment of the invention will be 

(42 and -43) divide the gray scale into three substantially described. Referring to FIG. 7, the first embodiment is an 

equal portions, as indicated in FIG. 5C. image-processing apparatus comprising an input image 

The four-by-four pattern matrix PM2, for example, has 4S memory 2, an address generator 4, a threshold processor 6, 

the value minus one hundred twenty-eight (-128) in all four a control unit 8 with an internal ring counter 8A, and an error 

cells in the upper left quadrant, the value one hundred diffusion processor 10. 

twenty-seven (127) in all four cells in the upper right The input image memory 2 receives and stores an exter- 

quadrant, the value forty-two (42) in all four cells of the nally supplied image signal with two hundred fifty-six gray 

lower left quadrant, and the value minus forty-three (-43) in 50 levels (from 0 to 255). The storage addresses correspond to 

all four cells of the lower right quadrant The four quadrants a Cartesian coordinate system, and will be expressed below 

of the other pattern matrices are likewise filled with these as (x\ y'). 

values - The address generator 4 receives Cartesian coordinates (x, 

Larger pattern matrices, such as thirty-two-by-thirty-two y) of pixels in the image from the control unit 8, and 

and sixty-four-by-sixty-four pattern matrices, can be simi- 55 generates corresponding memory addresses (x f , y'). These 

larly constructed and used for processing large images. memory addresses are supplied to the input image memory 

The reason for placing numbers with substantially equal 2, which supplies the corresponding pixel values to the 

absolute values in horizontally adjacent quadrants is to threshold processor 6 and error diffusion processor 10. The 

obtain an equal balance between positive and negative addresses are generated in a raster-scan sequence, starting in 

ofisets in each horizontal scanning line, similar to the equal 60 the top left corner of the image, following each horizontal 

positive and negative intervals of the sine functions used in scaiimng fine from left to right, proceeding downward from 
Fourier-s eries exp ansions. The reason for making the sum of » one horizontal scanning line to the next, and ending in the 

the four quadrants substantially equal to zero is to avoid bottom right comer. 

changing the overall darkness of the image. If the sum of the The error diffusion processor 10 comprises an adder 11, a 

four quadrants were much less than zero, the image would 65 comparator 12, a symbol generation unit 13, an image output 

be noticeably lightened; if the sum were much more than unit 14, a subtractor 15, an error memory 16, and a digital 

zero, the image would be darkened. filter 17. 
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The adder 11 adds each pixel value I received from the When the value in each quadrant of each table is also 

input image memory 2 to an error diffusion value received constant, as in the present case, the addressing process can 

from the digital filter 17, obtaining a modified pixel value I'. be further simplified as shown in FIG. 9. Each table stores 

The comparator 12 compares the modified pixel value I' onl y values. Table Tl is addressed by the least signifi- 

with a modified threshold Th' supplied by the threshold 5 cant b 1 lt , of each raonfaiule. ™>fc T2 is addressed I by the 

processor 6, outputs the maximum gray level CL- equal to S ^° nd ^ a f i ™* 13 . 15 

nzc\ e it • i » t . Vu, 7^ ♦ . tU addressed by the third-least significant bit of each coordi- 

255) if I' is equal to or greate fcan 1^ and outputs the flate ^ ^ [& &6drQs ^ £ fourth ., east si ificaQt Wt 

minimum gray level (1,^, equal to 0) otherwise. of eacfa Wben ^ ^ eight pixels £ „ ^ 

The symbol generation unit 13 converts the values (0, are read from the input image memory 2, for example, table 

<js 255) output by the comparator 12 to bit values or symbol T4 returns the same value eight times. 

values (0, 1). The values output from tables H, 12, 13, T4 will be 

The image output unit 14 temporarily stores the bit values denoted Kl, K2, K3, K4, respectively, 

generated by the symbol generation unit 13, and outputs Multiplier Ml multiplies the periodic function value Kl 

them in correspondence to image-coordinate information (x, 15 output from table Tl by the weight Wl output from the 

y). weighting table 31. Similarly, multiplier M2 multiplies the 

Tne subtractor 15 subtracts the value (1^ or l min ) output °*P* *2 of table T2 by weight W2, multiplier -M3 multi- 

by the comparator 12 from the modified valueT input P h f ^ output K3 of table 13 by ^weight .W3. and mulUplier 

. to the comparator 12 to obtain their difference, which is Ac M lmulUphes the output K4 of table T4 by weight W4 

A °\ quantization error of the pixel 20 r^JH 8 ^^^?^ ^T*^ 

U £K} pliers Ml, M2, M3, M4. The resulting sum of products Q is 

The error memory 16 temporarily stores the quantization similar to a Fourier series 

errors found by the subtractor 15 for the current horizontal Mul tiplier33 multiplies' the sum Q obtained by adder 32 

X scanning line and the preceding line. At any given mstant „ a bas£ coefficient „ ^ value ofaisa ^ 

^ the error memory 16 holds quantization error data for all „ ^ controluni , 8 and is approximately equal to 0.5 L the 

pucels to the left of the present pixel m the same scanning 25 embodiment. The output of multiplier 33 is theoffset 

line, and tor all pixels in tne scanning line immediately ATk ^ of 

multiplication bv a is to adiust the 

above. r A- : i J J 

range of onset values. 

The digital filter 17 selects four of the quantization error adder 34 adds me offsel va i ue AXh out t b 

values stored in the error memory 16, calculates a weighted me multipUer 53 to me basic threshold value (Th, equal to 

f™.°!L* if* leCtCd en0r ValUCS * 0Dd provides lbe result 10 128), which is the center value of the gray scale (halfway 

4i u„ 11 between 0 and 255). The output of adder 34 is the modified 



n/i 



the adder 11. 

The threshold processor 6 comprises a weighting table 31, threshold Th' supplied to the comparator 12 in the error 

further tables Tl, T2, 13, T4, multipliers Ml, M2, M3, M4, diffusion processor 10. 

an adder 32, another multipUer 33, and another adder 34. ^ xhe values of the weighting functions and the value of the 

The weighting table 31 stores the weighting functions base coefficient a should be selected according to the type of 

Wl, W2, W3, W4 described above. For each weighting output device (for example, a printer). The weighting func- 

function, the weighting table 31 stores weights from zero to tions do not have to employ sine or cosine curves as shown 

one as functions of gray levels from zero to two hundred in FIG. 6. Depending on the characteristics of the output < 

fifty-five, as depicted in FIG. 6. Alternatively, the sine curve ^ device, semicircular or triangular weighting functions can be' 

of weighting function W2 is stored, with supplementary used. The value of a can be decreased to a value less than 




logic for shifting this curve to the left or right to obtain the 0.5, such as 0.4, for example, to decrease the regularity of 

other weighting functions. The weighting table 31 receives the dot patterns in the output image, or increased to a value 

the gray level of a pixel from the input image memory 2, and greater than 0.5, such as 0.6, for example, to increase the 

outputs four corresponding weights Wl, W2, W3, W4. 4S regularity of the dot patterns. 

Tables Tl, 12, 13, T4 store four periodic functions built The embodiment shown in FIG. 7 comprises well-known 

from repetitions of the four pattern matrices PM1, PM2, computing hardware elements such as arithmetic and logic 

PM3, PM4 as described above. Each table is structured to processing elements and memory elements. Alternatively, 

1 receive a pair of Cartesian coordinates (x, y) and output the some of these elements can be replaced by software ele- 

f value of the corresponding periodic function at that coordi- 59 ments executed by a general-purpose processor such as a 

nate position. The table may simply store the function values microprocessor. 

at each coordinate position in the image plane. Alternatively, Next, the operation of the first embodiment will be 

the table may store one copy of the corresponding pattern described. 

matrix, supplemented by addressing logic for converting By way of introduction, first the conventional error dif- 

image coordinates (x, y) to positions in the pattern matrix. 55 fusion process will be explained through a simplified one- 

When the pattern matrices have sizes (numbers of pixels) dimensional example. FIG. 10 shows eight pixels having a 

that are powers of two, as in the present case, the coordinate constant gray level (I) slightly above a basic threshold value 

conversion can be carried out as shown in FIG. 8. Table Tl (Th). When these pixels are converted to bi- level pixels, the 

is addressed by the least significant bit of the x-coordinate, first pixel PI is made white (because I>Th), leaving a 

f and the least significant bit of the y-coordinate. Table T2 is 60 negative quantization error El. This error El is added to the 

addressed by two least significant bits of each coordinate. second pixel P2, which thus becomes black (because 

/ Table 13 is addressed by the three least significant bits of I+El<Th), leaving a positive quantization error E2. This 

each coordinate. Table T4 is addressed by the four least error E2 is added to the third pixel P3, and the above process 

significant bits of each coordinate. As the image is read from continues, producing white and black pixels as shown at the 

the input image memory 2, each position in each table is 65 bottom of FIG. 10. 

accessed repeatedly by means of the ring counter 8A in the In FIG. 10, the quantization error of each pixel is passed 

control unit 8. to the next pixel to the right. In two dimensions, the 
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quantization error is preferably diffused to several neigh- 
boring pixels. In the well-known Floyd-Steinberg error 
diffusion algorithm, illustrated in FIG. 11, seven sixteenths 
( 7 /i6) of the quantization error of pixel PI is added to pixel 
P1A, three sixteenths ( 3 /ie) to pixel P1B, five sixteenths (Vie) 
to pixel PIC, and one sixteenth (Vie) to pixel P1D. 

As a result, when an image is processed in raster-scan 
sequence, each pixel receives diffused quantization error 
from four previously processed pixels as shown in FIG. 12, 
in which an asterisk designates the pixel P2(x, y) being 
processed, E denotes quantization error, and A denotes a 
diffusion coefficient C 7 ' 16 - Vie. 5/10, or Me). 

The error diffusion process illustrated in FIGS. 11 and 12 
is carried out by the error diffusion processor 10 in the 
present embodiment. FIGS. 11 and 12 also illustrate the 
coordinate system used in the present embodiment, the 
positive x-direction extending to the right and the positive 
y^direction extending downward. 

The processing performed by the first embodiment will 
now be explained with reference to the flowchart in FIGS. 
13 and 14. 

First, the x and y coordinates are initialized to zero (step 
1), and the error memory 16 is initialized to all zero contents 
(step 2). These steps are performed by the control unit 8. 

The following steps form a loop that is repeated once for 
each pixel in the image. 

First, the digital filter 17 calculates the quantization error 
E(x, y) to be added to the pixel at position (x, y), (step 3). 
As illustrated in FIG. 12, E(x, y) is calculated by the 
following weighted sum. 

EQt, y^A(x^y-iyE(x-^y-l) *AQc, y-l)E(x, jy-l) +A(*+l f 
y-l)E(x+l, y-1) 4*4(*-l, y)-E(x-l, y) 

Next, the adder 11 adds E(x, y) to the pixel value I(x, y) 
to obtain a modified pixel value I'(x, y) (step 4). 

Next, the values Kj(x, y) of the four periodic functions 
(j=l, 2, 3, 4) at coordinates (x, y) are obtained from tables 
Tl, 12, T3, T4 (step 5). 

These values Kj(x, y) are multiplied by corresponding 
weights Wj(I(x, y)) obtained from the weighting table 31, 
and added by adder 32 to obtain the sum Q. For brevity, Wj 
will be written instead of Wj(I(x, y)) (j=l, 2, 3, 4). 

Q-K1(& y) TKl+K2(x; y) Wl^K3(x, y)Wy+K4& y)-WA 

Multiplier 33 multiplies this sum Q by the base coefficient 
a to obtain the offset ATh(x, y), and adder 34 adds the offset 
ATh(x, y) to the basic threshold Th to obtain the modified 
threshold Th'(x, y) (step 6). 

Next, the comparator 12 compares the modified threshold 
Th'(x, y) with the modified pixel value I'(x, y) (step 7). 

If the modified threshold Th'(x, y) is less than the modi- 
fied pixel value I'(x, y), then the symbol generation unit 13 
outputs the symbol "1" (step 8), and the subtractor 15 
calculates the quantization error E(x, y) as the difference 
I'(X, yHm** (step 9). If the modified threshold Th'(x, y) is 
equal to or greater than the modified pixel value I'(x, y), then 
the symbol generation unit 13 outputs the symbol "0" (step 
10), and the subtractor 15 calculates the quantization error 
E(x, y) as the difference I'(x, y)-\ min (step 11). The quanti- 
zation error E(x, y) calculated in step 9 or step 11 is stored 
in the error memory 16. 

Incidentally, the same notation E(x, y) is used in FIG. 13 
to denote both the quantization error diffused into position 
(x, y) in steps 3 and 4, and the quantization error diffused out 
from position (x, y) in steps 9 and 11. 
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Referring to FIG. 14, following step 9 or step 11, the 
x-coordinate (x) of the pixel that has just been processed is 
compared with the maximum x-coordinate (x^) of the 
image (step 12). If x is less than x^, then x is incremented 

5 (step 13), and the process returns to step 3 to process the next 
pixel in the same horizontal scanning line. If x is not less 
than x maxy then the y-coordinate (y) of the pixel is compared 
with the maximum y-coordinate (y mfl J of the image (step 
14). If y is less than y,^ then x is set to zero and y is 

1Q incremented (step 15), and the process returns to step 3 to 
process the first pixel in the next horizontal scanning line. If 
y is not less than y ma3fi the process ends. 

Steps 12 to 15 form a decision procedure that determines 
whether the entire image has been processed, proceeds to the 
next pixel if the entire image has not been processed, and 

15 otherwise terminates the processing, to await a command to 
process the next image. 

The process described above generates a bi-level image in 
which the spatial frequencies of an area are constrained by 
establishing regular, hierarchical frameworks having posi- 

20 tions where dots are encouraged to appear, and positions 
where dots are encouraged not to appear, responsive to the 
gray level represented by the pseudo-tone pattern of the area. 
Although a certain degree of randomness still appears in the 
pseudo-tone patterns, the randomness is constrained; the 

25 output dots have a strong tendency to appear at or near 
preferred positions, and the dot patterns show large-scale 
regularity. The dot noise produced by conventional error 
diffusion is thereby reduced. In a printer, dot positions can 
be controlled according to the characteristics of the printer, 

30 to obtain an output image of high quality. 

The pattern matrices from which the periodic functions 
are generated push pixels toward white in the second and 
fourth quadrants, and toward black in the first and third 
quadrants. The output dot patterns are thereby given a 

35 desirable checkerboard appearance. The quadrants of the 
pattern matrices expand hierarchically from the first periodic 
function to the fourth periodic function, and the dot patterns 
show similar hierarchical regularity. 
Although the first embodiment uses an error diffusion 

40 processor 10, this can be replaced by a direct dither pro- 
cessing unit of the type that employs a dither matrix. The 
threshold processor 6 then offsets each of the threshold 
values appearing in the dither matrix. That is, a matrix of 
fixed thresholds is replaced with a matrix in which the 

45 thresholds are expressed as Fourier series. 

Next, a second embodiment will be described. 
Referring to FIG. 15, the second embodiment inserts an 
accumulator 19 between the subtractor 15 and error memory 
26, a coordinate converter 20 between the control unit 28 

50 and address generator 4, and a divider 21 between the digital 
fitter 17 and adder 11. 

The second embodiment performs error diffusion on a 
block-by-block basis. Referring to FIG. 16, the input image 
is divided into rectangular m-by-n blocks, where m and n are 

55 positive integers. The entire image is covered by these 
rectangular blocks. 

The control unit 28 controls the image processing so that, 
instead of processing one horizontal scanning line at a time, 
the error diffusion processor 30 processes one block at a 

60 time. First the block in the upper-left corner is processed, the 
pixels within this block being selected in raster-scan 
sequence, from left to right and from top to bottom. Then the 
next block to the right is similarly processed. When the right 
edge of the image is reached, the leftmost block in the next 

65 row of blocks is processed 

In the description below, blocks are identified by block 
coordinates (xb, yb) and pixels within a block are identified 
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by an index (i). In the processing of the first block, for 
example, xb and yb are zero, and i varies from one (1) to a 
limit L equal to the number of pixels per block (L=mn). 

Id each block, the accumulator 19 calculates the cumu- 
lative sum of the quantization errors supplied by the sub- 5 
tractor 15. At the end of the block, the cumulative sum is 
stored in the error memory 26, and the accumulator 19 is 
reset by the control unit 28 in preparation for the next block. 

The error memory 26 thus stores one error value per 
block, instead of one error value per pixel. At any given 10 
moment, the error memory 26 contains error values for the 
blocks to the left of the block now being processed in the 
same row, and for all blocks in the row just above. 

The digital filter 17 performs the same type of calculation 
as in the first embodiment, using block coordinates and 15 
block error values instead of pixel coordinates and pixel 
error values. As shown in FIG. 17, each block receives error 
diffusions from the three neighboring blocks in the row 
above and the block immediately to the left in the same row. 
The diffusion coefficients Ab are, for example, the same as 20 
the coefficients A in FIG. 12. 

The coordinate converter 20 converts the block coordi- 
nates (xb, yb) and pixel index value (i) output by the control 
unit 28 to pixel coordinates (x, y). 

The divider 21 divides the output of the digital filter 17 by 25 
the number of pixels per block (L), thereby obtaining the 
per-pixel error diffusion. If L is a power of two, the divider 
21 can be configured as a bit shifter. 

The threshold processor 6 performs the same processing 
as in the first embodiment, obtaining a separate threshold 30 
offset ATh(x, y) for each pixel. The order in which the pixels 
are processed differs from the first embodiment, but the 
offset values are the same. 

Next, the operation of the second embodiment will be 
described with reference to the flowchart in FIGS. 18 and 19. 35 

To begin the processing of an image, the block coordi- 
nates (xb, yb) are cleared to zero (step 101) and the error 
memory 26 is initialized (step 102). 

Next, the digital filter 17 calculates the total error Eb(xb, 
yb) diffused into block (xb, yb) as follows (step 103), using ^ 
the block error values (Eb) of neighboring blocks stored in 
the error memory 26, and diffusion coefficients (Ab) as 
shown in FIG. 17. 

Eb{xb> yb)«Ab(xb-l, yb-1)Eb{xb-\, yb-l) +Ab{xb, yb-iyEb{xb, 
yb-1) +Ab(xb+1 , yb-l)Eb(xb+l , yb-l) +Ab(xb-1 , yb)Eb(xb- 45 
1.3*) 

The pixel index (i) is now initialized to one, indicating the 
pixel in the top left comer of block (xb, yb) (step 104). 

The adder 11 modifies the pixel value I(xb, yb, i) by 
adding the quotient of Eb(xb, yb) divided by the number of 50 
pixels (L) in the block, obtaining a modified pixel value 
I*(xb, yb, i) (step 105). 

Next, the values Kl, K2, K3, K4 of the four periodic 
functions at coordinates (x, y) are obtained from tables Tl, 
T2, T3, T4 (step 106). The coordinates (x, y) are generated 55 
by the coordinate converter 20, and correspond to the block 
coordinates and index (xb, yb, i) output by the control unit 
28. 

Next, the periodic function values Kl, K2, K3, K4 are 
multiplied by corresponding weights Wl, W2, W3, W4 60 
obtained from the weighting table 31, and the products are 
added by adder 32 to obtain a sum Q, as in the first 
embodiment. Multiplier 33 multiplies the sum Q by the base 
coefficient a to obtain the offset ATh(xb, yb, i) (step 107). 

Adder 34 adds the offset ATh(xb, yb, i) to the basic 65 
threshold Th to obtain the modified threshold Th'(xb, yb, i) 
(step 108). 
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Referring to FIG. 19, the comparator 12 compares the 
modified threshold Th'(xb, yb, i) with the modified pixel 
value T(xb, yb, i) (step 109). 

If the modified threshold Th'(xb, yb, i) is less than the 
modified pixel value I'(xb, yb, i), then the symbol generation 
unit 13 outputs the symbol "1" (step 110), and the subtracter 
15 calculates the quantization error E(xb, yb, i) of the pixel 
just processed as the difference I'(xb, yb, ij-l^ (step 111). 
If the modified threshold Th'(xb, yb, i) is equal to or greater 
than the modified pixel value F(xb, yb, i), then the symbol 
generation unit 13 outputs the symbol "(T (step 112), and the 
subtracter 15 calculates the quantization error E(xb, yb, i) as 
the difference I'(xb, yb, $-\ min (step 113). 

Next, the pixel index (i) is incremented (step 114), and 
compared with the number (L) of pixels per block (step 115). 
If the incremented pixel index is equal to or less than L, the 
process returns to step 105 to process the next pixel in the 
block. 

As steps 105 to 115 are repeated for each pixel in the 
block, the accumulator 19 calculates a cumulative sum 
2E(xb, yb, i) of the quantization errors of the pixels. When 
the incremented pixel index exceeds L, this cumulative sum 
is stored in the error memory 26 as the block error Eb(xb, 
yb) (step 116). The same notation Eb(xb, yb) is used in 
FIGS. 18 and 19 to denote both the quantization error 
diffused into block (xb, yb) in step 103, and the quantization 
error diffused out from block (xb, yb) in step 116. 

Next, the x-coordinate (xb) of the block that has just been 
processed is compared with the maximum block 
x-coordinate (xb^J in the image (step 117). If xb is less 
than xb majcy then xb is incremented (step 118), and the 
process returns to step 103 to process the next block in the 
same horizontal row. If xb is not less than xb^^ then the 
y-coordinate (yb) of the block is compared with the maxi- 
mum block y-coordinate (yb^^ in the image (step 119). If 
yb is less than yb^^, then xb is set to zero and yb is 
incremented (step 120), and the process returns to step 103 
to process the first block in the next row. If yb is not less than 
ybmax, the process ends. 

Steps 117 to 120 cause all blocks to be processed in turn, 
until the entire image has been processed. 

The second embodiment has the same effect as the first 
embodiment in imposing a desirable regularity on the dot 
patterns used to represent each gray level, thus reducing dot 
noise in the converted image. In addition, the amount of 
randomness appearing in the output dot patterns decreases as 
the block size increases, providing a method of controlling 
the amount of randomness. 

A further advantage of the second embodiment is that the 
digital filter 17 has to operate only once per block, instead 
of once per pixel, reducing the filtering- processing load by 
a factor of L. Accesses to the error memory 26 are also 
reduced by a factor of L. The error diffusion processor 30 in 
the second embodiment thus operates more efficiently than 
the error diffusion processor 10 in the first embodiment. 

In a variation of the second embodiment, the quantization 
error diffused into a block is distributed in a weighted 
manner among the pixels in the block, according to the pixel 
position, for example, instead of being distributed to all 
pixels equally. 

Next, a third embodiment will be described. The third 
embodiment reduces the gray scale of the input image to 
four output levels instead of two. 

The four output levels are denoted gO (black), gl (dark 
gray), g2 (light gray), and g3 (white) in FIG. 20. There are 
now three basic thresholds: Thl midway between gO and gl, 
Th2 midway between gl and g2, and Th3 midway between 
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g2 and g3. Each of these basic thresholds is modified by the in each color plane are controlled *so that dots of different 

offset ATh calculated by the threshold processor 6 for each colors do not overlap and clumping is avoided Specifically, 

pixel. The threshold processor 6 operates in the same way as the phase of the periodic functions is varied from one color 

in the first or second embodiment, except that the value of plane to another, by providing the control signal (j) as a third 

the base coefficient a is preferably smaller than in those 5 m P ul t0 Tl, TO, T3, T4, by means of signal lines not 

embodiments, and the second adder 34 adds ATh to each of shown m mG - 21 

the basic thresholds Thl, Th2, Th3 separately, generating Conversely, when the C, M, and Y values are large, the dot 

three modified thresholds Thl', Th2', Th3'. patterns of different colors are controlled so as to overlap, by 

The error diffusion processor 10 in the third embodiment al £T g P^es of the periodic : functions in each plane, 

is the same as the error diffusion processor 10 in the first 10 which results in better rendition of pseudo-tone patterns in 

embodiment, with the following modifications. The com- SZbe differences described above, in the pro- 

^ ii Tftf^ ™ ^ each color plane, the threshold processor 36 P and 

adder 11 with the modified thresholds Thl', Th2\ TWw error ^ 10 te ^ m ^ e fifSt embodi . 

turn to select one of the four gray levels gO, gl, g2, g3. The mcnt 

symbol generation unit 13 outputs symbols "0" "1," «2," "3" is fourth embodiment provides the same advantage as 

corresponding to these gray levels. The subtracter 15 sub- the preceding embodiments in reducing dot noise, with the 

tracts the selected gray level from the modified pixel value additional advantage of enabling the dot patterns in different 

to obtain the quantization error. color planes to be controlled so that the dot patterns are 

In a variation of the third embodiment, the block error overlaid on one another in desirable ways, 

diffusion processor 10 of the second embodiment is 20 In a variation of the fourth embodiment, block error 

employed, with similar modifications. diffusion is employed as in the second embodiment. 

In the third embodiment, a gray level higher than thresh- In another variation, the number of output levels in each 

old Th3 but lower than level g3, for example, is expressed color plane is increased to a number greater than two, as in 

as a pattern of dots with levels g2 and g3, as indicated at the me third embodiment. 

bottom of FIG. 20. The same benefits as in the first and 25 ^ described above, the first embodiment performs a 

second embodiments are obtained in terms of increased constrained error diffusion process that produces a bi-level 

regularity of the dot patterns expressing each gray level in ima e c ^ r^^ccd dot noise. The second embodiment 

the input image. performs a constrained block error diffusion process that 

Other variations of the third embodiment can be obtained ako produces a bi-level image with reduced dot noise. The 

by changing the number of output levels from four to three, 30 embodiment extends these processes to output images 

five, or any larger number N. w i* n more tnan two levels. The fourth embodiment extends 

Next, a fourth embodiment will be described. The fourth processes to color images, permitting control over the 

embodiment processes a color image. wav in which dot patterns in different color planes are 

Referring to FIG. 21, the input image in the fourth overlaid. All of these embodiments enhance the quality of 

embodiment comprises cyan, magenta, and yellow color 35 images printed by dot printing devices such as ink-jet 

planes. Three input image memories 2C, 2M, 2Y are pro- printers. 

vided for storing the pixel values in these three color planes. h * anticipated that the invention can also enhance the 
The addresses (x\ y) generated by the address generator 4 appearance of images displayed on various types of dot- 
are supplied to all three input image memories 2C, 2M, 2Y matrix display devices. 

simultaneously. A switch 22 selects the output from one 40 Those skilled in the art will recognize that further v aria- 
input image memory, under control of a signal (j) from the tions are possible within the scope claimed below, 
control unit 8. First, input image memory 2C is selected and What is claimed is: 

the cyan plane of the entire image is processed. Next, input 1 - A memod of processing input image data to generate 

image memory 2M is selected and the magenta plane of the pseudo-tone patterns, comprising the steps of: 

entire image is processed. Finally, input image memory 2Y 45 ( a ) modifying at least one predetermined threshold by 

is s elected and the yellow plane of the entire image is combining said predetermined threshold simulta- 

processed. neously with at least two periodic functions having 

The threshold processor 36 in the fourth embodiment has different periods, said periodic functions being func- 

an additional divider 35 that receives the pixel values C, M, ^ons of position in said input image data; 

Y output from the three input image memories 2C, 2M, 2Y, 50 (b) comparing said input image data with the modified 

and the pixel value Z selected by the switch 22 (Z-C, M, or threshold; and 

Y). The divider 35 calculates the following quantity , which (c) generating output image data with N levels, using 

is supplied to the weighting table 31. results obtained in said step (b), levels present in the 

j . z/ ( C+M+y) m P ut image data but not present among said N levels 

55 being simulated by pseudo-tone patterns in said output 

When the output image is printed by an ink-jet printer, and image data with N levels, N being an integer greater 

an arbitrary color is simulated by overlaid patterns of cyan, than one. 

magenta, and yellow dots, the quotient Z/(C+M+Y) indi- 2. The method of claim 1, wherein said input image data 

cates the share of the selected color plane in the total amount represent a plurality of horizontal scanning lines, and said 

of ink. The weights Wl, W2, W3, W4 used in each color 60 periodic functions are one-dimensional, following said hori- 

plane are thereby controlled according to the relative quan- zontal scanning lines. 

tity of ink printed in each of the three primary colors. 3. The method of claim 1, wherein said periodic functions 

Accordingly, the dot pattern of each one of the three primary are two-dimensional. 

colors is controlled in a way that takes the space occupied by 4. The method of claim 1, wherein said step (c) comprises 

the other two primary colors into account 65 error diffusion. 

In addition, when the values C, M, Y obtained from the 5. The method of claim 4, wherein said step (c) further 

input image memories 2C, 2M, 2 Y are small, the dot patterns comprises the steps of: 
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(d) dividing said input image data into blocks; 

(e) calculating a total quantization error in each of said 
blocks; and 

(f) distributing weighted portions of said total quantiza- 
tion error to neighboring blocks. 

6. The method of claim 1, wherein N is equal to two. 

7. The method of claim 1, wherein N is greater than two. 

8. The method of claim 1, wherein said input image data 
are color image data comprising a plurality of color planes, 
and said steps (a), (b), (c) are applied separately to each of 
said color planes. 

9. A method of processing input image data to generate 
pseudo-tone patterns, comprising the steps of: 

(a) modifying at least one predetermined threshold 
according to at least two periodic functions having 
different periods, said periodic functions being func- 
tions of position in said input image data by multiply- 
ing each of said periodic functions by a weighting 
function to obtain products, using different weighting 
functions for different periodic functions, said weight- 
ing functions being functions of the levels in said input 
image data, and adding the products to obtain a sum; 

(b) comparing said input image data with the modified 
threshold; and 

(c) generating output image data with N levels, using 
results obtained in said step (b), levels present in the 
input image data but not present among said N levels 
being simulated by pseudo-tone patterns in said output 
image data with N levels, N being an integer greater 
than one. 

10. The method of claim 9, wherein said step (a) further 
comprises the steps of: 

(i) multiplying said sum by a constant coefficient to obtain 

an offset; and 
(j) adding said offset to said predetermined threshold. 

11. The method of claim 9, wherein each said weighting 
function has non-zero values over a predetermined range of 
levels, and is zero outside said predetermined range. 

12. The method of claim 11, wherein said non-zero values 
follow a sine curve within said predetermined range. 

13. The method of claim 12, wherein as the periods of said 
periodic functions decrease, the ranges of non-zero values of 
the weighting functions of said periodic functions shift 
toward higher levels. 

14. The method of claim 13, wherein the ranges of 
non-zero values of said weighting functions mutually over- 
lap. 

15. An image-processing apparatus for processing input 
image data to generate output image data with pseudo-tone 
patterns, comprising: 

a threshold processor modifying at least one predeter- 
mined threshold by combining said predetermined 
threshold simultaneously with at least two periodic 
functions having different periods, said periodic func- 
tions being functions of position in said input image 
data; and 

an error diffusion processor coupled to said threshold 
processor, comparing said input image data with the 
modified threshold, and using comparison results thus 
obtained to generate output image data with N levels, 
in which levels present in the input image data but not 
present among said N levels are simulated by pseudo- 
tone patterns in said output image data with N levels, 
N being an integer greater than one. 



15 
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16. The image-processing apparatus of claim 15, wherein 
said input image data represent a plurality of horizontal 
scanning lines, and said periodic functions are one- 
dimensional, following said horizontal scanning fines. 

17. The image-processing apparatus of claim 15, wherein 
said periodic functions are two-dimensional. 

18. The image-processing apparatus of claim 15, wherein 
said error diffusion processor performs an error diffusion 
process. 

19. The image-processing apparatus of claim 18, wherein 
said error diffusion processor divides said input image data 
into blocks, calculates a total quantization error in each of 
said blocks, and distributes weighted portions of said total 
quantization error to neighboring blocks. 

20. The image-processing apparatus of claim 15, wherein 
N is equal to two. 

21. The image-processing apparatus of claim 15, wherein 
N is greater than two. 

22. The image-processing apparatus of claim 15, wherein 
said input image data are color image data comprising a 
plurality of color planes, and said threshold processor and 
said error diffusion processor operate on each of said color 
planes separately. 

23. An image-processing apparatus for processing input 
image data to generate output image data with pseudo-tone 
patterns, comprising: 

a threshold processor modifying at least one predeter- 
mined threshold according to at least two periodic 
functions having different period, said periodic func- 
tions being functions of position in said input image 
data; 

an error diffusion processor coupled to said threshold 
processor, comparing said input image data with the 
modified threshold, and using comparison results thus 
obtained to generate output image data with N levels, 
in which levels present in the input image data but not 
present among said N levels are simulated by pseudo- 
tone patterns in said output image data with N levels, 
N being an integer greater than one; 

wherein said threshold processor multiplies each of said 
periodic functions by a weighting function, using dif- 
ferent weighting functions for different periodic 
functions, said weighting functions being functions of 
the levels in said input image data, and adds resulting 
products to obtain a sum. 

24. The image-processing apparatus of claim 23, wherein 
said threshold processor multiplies said sum by a constant 
coefficient to obtain an offset, and adds said offset to said 
predetermined threshold. 

25. Th e image-processing apparatus of claim 23, wherein 
each said weighting function has non-zero values over a 
prea^termined range of levels, and is zero outside said 
predetermined range. 

26. The image-processing apparatus of claim 25, wherein 
said non-zero values follow a sine curve within said prede- 
termined range. 

27. The image-processing apparatus of claim 26, wherein 
as the periods of said periodic functions decrease, the ranges 
of non-zero values of the weighting functions of said peri- 
odic functions shift toward higher levels. 

28. The image-processing apparatus of claim 27, wherein 
the ranges of non-zero values of said weighting functions 
mutually overlap. 
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